Number variance for hierarchical random walks 
and related fluctuations* 

Tomasz Bojdecki 1 Luis G. Gorostiza^ Anna Talarczyk 1 

tobojd@mimuw.edu.pl lgorosti@math.cinvestav.mx annatal@mimuw.edu.pl 



Abstract 

We study an infinite system of independent symmetric random walks on a hierarchical 
group, in particular, the c-random walks. Such walks are used, e.g., in population genetics. 
The number variance problem consists in investigating if the variance of the number of 
"particles" N n (L) lying in the ball of radius L at a given time n remains bounded, or 
even better, converges to a finite limit, as L — > oo. We give a necessary and sufficient 
condition and discuss its relationship to transience/recurrence property of the walk. Next 
we consider normalized fluctuations of N n (L) around the mean as n — > oo and L is 
increased in an appropriate way. We prove convergence of finite dimensional distributions 
to a Gaussian process whose properties are discussed. As the c-random walks mimic 
symmetric stable processes on R, we compare our results to those obtained by Hambly and 
Jones (2007,2009), where the number variance problem for an infinite system of symmetric 
stable processes on R was studied. Since the hierarchical group is an ultrametric space, 
corresponding results for symmetric stable processes and hierarchical random walks may 
be analogous or quite different, as has been observed in other contexts. An example of a 
difference in the present context is that for the stable processes a fluctuation limit process 
is a centered Gaussian process which is not Markovian and has long range dependent 
stationary increments, but the counterpart for hierarchical random walks is Markovian, 
and in a special case it has independent increments. 
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I 



1 Introduction 



Random walks on hierarchical groups, called hierarchical random walks, have arisen in several 
applications. The hierarchical random walks we will consider in this paper are of the type 
introduced by Spitzer [21] (p. 93) in a special case, and more generally by Sawyer and Felsenstein 
[20J in the context of population genetics. We will focus on "number variance" properties and 
related fluctuations for these random walks. This question has been investigated by Hambly 
and Jones P3J [H] for cv-stable processes on R. 
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t Centro de Investigation y de Estudios Avanzados, A. P. 14-740, Mexico 07000 D.F., Mexico. 
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For an integer M > 2, the hierarchical group of order M is defined by 

= {X = (Xi, X 2 , ■■■): Xi e {0,1, ... ,M - 1}, HiXi < oo}, 

with addition componentwise mod M. It is a countable Abelian group which is also described as 
the direct sum of a countable number of copies of the cyclic group of order M. The hierarchical 
distance | • | on is defined by 

I | _ J if x — y, 

|X 2/1 \ max{i : 7^ y^} if x 7^ y. 

It is translation-invariant, and satisfies the strong (non- Archimedean) triangle inequality 

|x — y| < max{|x — z\, \z — y\} for any x,y,z. 

This means that (CIm, | • |) is an ultrametric space. In such a space two balls are either disjoint 
or one is contained in the other. This property is the cause of significant differences from 
mathematical models based on Euclidean spaces R d or Euclidean lattices Z d , and it is interesting 
to study analogies and differences between "corresponding" models defined on Q,m and on R d 
or Z d , specially because hierarchical models may reveal behaviors that are not observed in 
Euclidean models. A picture of (£Im, I • I) is the set of leaves at the top of an infinite tree where 
each inner node at each level j > 1 has one predecessor at level j + 1 and M successors at level 
j — 1. The distance between two individuals (leaves) at level j = is the depth in the tree to 
their most recent common ancestor, and it measures the degree of relatedness between the two 
individuals. 

A random walk £ = (£„) n >o on Om starting at is defined by £0 = 0, £ n = p\ + - • -+p n , n > 1, 
where pi, p 2 , . . . are independent copies of p, which is a random element of VLm with distribution 
of the form 

p{p = v) = m^W^T)' p (" = > = °. < L1 > 

where (7^)^=1,2,... is a probability law on {1,2,...}. That is, the jumps of £ n are taken by first 
choosing distance j with probability r-j, and then choosing a point with uniform probability 
among those at distance j from the previous position of the walk (note that M j ~ 1 (M — 1) is the 
number of points at distance j from a given point of ^m)- These random walks are the most 
general symmetric random walks on Qm- We assume that (rj) is not restricted to a bounded 
set, so that the walk can reach arbitrarily large distances (by ultrametricity, it is not possible 
to go far with small steps). We refer to a hierarchical random walk determined by (r-j) as rj-rw. 
A particular r^-rw is the c-rw, where c is a constant such that < c < M, and 

The c-rw mimics the behavior of the standard a-stable process on M d in the sense that both have 
analogous recurrence/transience behaviors. The analogies are given in terms of their degrees 7. 
The Appendix contains information on 7. In particular, 7 < corresponds to recurrence and 
7 > to transience. The a-stable process on M. d has degree 

7 = -"l, (1-3) 
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and the c-rw on Vt M has degree 

7 log(M/c)' 1 ' J 

By relating the degrees (11.31) and fll.4f) we will compare our results with those of [131 E] • Other 
jj-rw's we will consider are defined in the Appendix and used for examples. 

Hierarchical structures are found in the physical, biological and social sciences due to the 
multiscale organization of many natural objects (see e.g. [21 [19]). Some of the main applications 
of hierarchical systems are found in statistical physics; a partial list of references is contained 
in [7j. We have already mentioned [20] in population genetics. Some other references where 
stochastic models based on hierarchical structures have been studied are [101 (121 EH] (inter- 
acting diffusions), [3 [8] (branching systems), [1] (contact processes), [31 [H [18] (percolation), 
[T5l [16] (search algorithms) . 

The "number variance" problem studied [131 [H] f° r a-stable processes on K has a motivation 
from physics related to so-called determinantal processes (see the Introduction of [13])- The 
model is as follows in simplified form. An independent copy of the standard a-stable process, 
representing the motion of a particle, starts at time from each point Uj = j — e,j G Z, where 
e is a uniformly distributed random variable on [0, 1] (the e-displacement provides spatial 
homogeneity). Let Nt[0, L] denote the number of particles that at time t > lie in the interval 
[0, L], and consider its variance, VarA^[0, L]. The "number variance" problem refers to the 
behavior of VaiN t [0, L] for fixed t as L — > oo. It is shown that for a < 1 (transient case) the 
variance tends to oo, and for 1 < a < 2 (recurrent case) it has a finite limit up to an additive 
fluctuation bounded by 1. The latter behavior is called saturation. For a = 1 the process is 
recurrent and the variance tends to infinity. The rescaled fluctuation process defined by 

w = »y^ | s£0i (15) 

is discussed, and a continuous interpolation of this process is shown to converge weakly on 
a space of continuous functions, as t — > oo, to a centered Gaussian process with stationary 
increments, which is not Markovian and has long range dependence. Hence this process has a 
resemblance to fractional Brownian motion. 

The purpose of the present paper is to prove counterparts to the results of [131 EH] f° r a 
system of independent hierarchical random walks starting from each point of Qm- We will do 
this for a general class of hierarchical random walks, and we will also consider a model where 
at the initial time there is a Poisson number of particles at each site. As with other systems 
of hierarchical random walks, one of our objectives is to exhibit analogies and differences with 
the results for the Euclidean model. 

Now we give a summary of our results. 

We consider a system of independent r 3 -rw's (particles) starting from each point of Qm, 
and denote by N n (L) the number of particles that at step n lie in the ball {x G Qm '■ \%\ < L}, 
L G Z + (in [T3"l [H] the interval [0, L] was considered, but essentially nothing changes if it is 
replaced by [— L, L]). For general rj we give a necessary and sufficient condition for boundedness 
of Var N n (L) and for existence of its limit as L — > oo. It turns out that 

lim Var iV^L) = 2n— - — lim M L r L 
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(Theorem 2.2). In particular, for the c-rw this limit exists, and it is finite if and only if 
c < 1, i.e., if the walk is recurrent. For the non-critical cases, c < 1 and c > 1, this result 
corresponds to the properties of the number variance for a-stable processes mentioned above, 
but the correspondence breaks down for the critical case, c = 1: the process is recurrent and 
the variance has a finite limit. For general rj, boundedness of VariV n (L) and recurrence are 
not equivalent (see Remark 2.4(c)). For the initial Poisson system the situation is simpler, 
VariV n (L) always tends to infinity. It can be shown that the same thing happens for the 
a-stable process. 

Next we investigate the fluctuations of N n (L) around the mean as n — > oo, when L = L(n) 
is increased in an appropriate way. Analogously as in fll.5p . where a "time" parameter s was 
introduced, we introduce a "time" parameter t, taking L(n) + R(t), where R(t) is a suitable 
non-decreasing function of t, and we consider a corresponding fluctuation process. We assume 
that linij^oo rj+i/rj = a > 0. It turns out that the cases a < 1 and a = 1 are significantly 
different (see Theorem 2.8). For a < 1, the most appropriate L(n) has logarithmic growth 
and the norming is of the order n 9 ^ 2 , where 9 = log M/ log-. In general the fluctuation 
processes may not converge, but convergent subsequences can be chosen (in the sense of finite 
dimensional distributions), and a family of processes indexed by a parameter k G [1, -] is 
obtained. These processes have stepwise trajectories, they are Gaussian with covariance of the 
form f(sVt)g K (s At), in particular they are Markov. If a = 1, the fluctuation process converges 
to a Gaussian process with independent increments. In contrast, for the a-stable process on R 
the fluctuation limit is a non-Markov process with stationary increments [131 EH]- 

We study some further properties of our limit processes and discuss a probabilistic inter- 
pretation of the parameter 9, which is closely related to the degree 7. 

Finally, for the Poisson model we show that with general rj the fluctuation limit exists and 
it is, up to a constant, the same process as that obtained in the previous case for a = 1. 

Due to the ultrametric structure of Om, the calculations involved in the proofs are quite 
different from those in [T3~| [14] . 



2 Results 

Fix an integer M > 2 and let Qm be the corresponding hierarchical group as defined in Intro- 
duction. We consider a system of independent rj-rw's on Qm starting from each point u G Qm- 
Recall that N n (L) is the number of particles lying in the ball Bl = {iG Qm '■ M < L} at step 
n, i.e., 

N n (L)= + (2-1) 

where {(, u } u en M are independent copies of r^-rw's starting at 0. 
Proposition 2.1 

(a) EN n (L) = M L , n, L = 0,1,2, (2.2) 
(6) VariV n (L) ~ 2nM L P(\p\ > L) as L ->■ 00, (2.3) 
(c) P(\p\ < L) 2n ~ 1 nP{\p\ > L)M L < VariV ra (L) < 2M L , n, L — 0,1,2, ... , (2.4) 

(~ means that the quotient of both sides tends to 1). 
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The number variance problem is solved in the following theorem. 



Theorem 2.2 For any n — 1, 2, ... , limsup^^ VariV ra (L) < oo if and only if 

lim supi-j.oo M l tl < oo. Moreover, lirn^oo VariV n (L) exists (finite or not) if and only if 

lim^oo M L r L exists, and in this case 

lim VarAUL) = 2n 1 - lim M L r L . (2.5) 

This theorem applied to the examples mentioned in the Appendix gives the following results. 
Example 2.3 (a) c-rw: 



(b) jP-rw: 



lim VariV„(L) 

L—too 



lim VarAUL) = < ., V, ' 



i 


f c<l, 




2n i 


f c=l, 


(2.6) 


oo i 


f c> 1. 





^/or (3 = it is obviously a c-rw with c = 1 ) 
(c)r j = Dj-f > ,0>l: 

lim VariV n (L) = oo. 



Remark 2.4 (a) Let us compare the result (12. 6p for the c-rw with the solution of the number 
variance problem obtained in [13,14] for the a-stable process. With c < 1 (recurrent case) we 
obtain "true" limits, whereas for the a-stable process with a > 1 the variances are bounded 
but do not converge due to an oscillating term (see [14]). On the other hand, with c < 1 the 
limit is trivial (zero), whereas for the a-stable process with a > 1 the variances are bounded 
away from 0. The only non-trivial limit for the c-rw is obtained in the critical case c = 1, and 
it has no counterpart for the a-stable process, since for a = 1 the variances tend to infinity. 

(b) For the c-rw, finiteness of the variance limit is equivalent to recurrence, whereas for the 
a-stable process this equivalence breaks down in the critical case a = 1. 

(c) Example 2.3(b) shows that in general finiteness of the variance limit is not equivalent to 
recurrence (see the Appendix). 

We now give a result for the Poisson case. 

Proposition 2.5 Assume that initially at each site there is a Poisson number of particles, and 
these numbers are i.i.d. Then for any system of rj-rw's, lim^oo Var iV n (L) = oo. 

Remark 2.6 (a) It can be shown that an analogous result holds for a-stable processes, 
(b) We will come back to the Poisson system later on. Poisson systems seem to be the most 
natural as random initial configurations. However, it can be shown that for each initial config- 
uration determined by i.i.d. random variables {v x }xen M ( l/ x particles at site x) which are truly 
random, i.e., with Var v x > 0, we have lim^oo VariV n (L) = oo. 
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So far we have considered N n (L) for n fixed. Now we will vary both n and L, more 
precisely, we want to investigate the normalized fluctuations of N n (L) as n — > oo, simultaneously 
increasing L in an appropriate way. 

We make the following assumption on the random walk: 

r j+i < r j for j > some jo, and lim = a > 0. (2.7) 

This assumption is satisfied for all our examples (see Example 2.3 and the Appendix). 
By Proposition 2.1(c), it is natural to take 

L{n) = sup |l £ Z + : h(L) — ~ J > ( 2 - 8 ) 

where 

oo 

h(L) = P(\p\ > L) = r v (2-9) 

j=L+l 

so that, by (12 .4p . VaxN n (L(n)) has the same rate of increase as M L{ - n \ Hence VM L ^ is the 
natural normalization for the fluctuations of N n (L(n)). 
L(n) has the following properties. 



Lemma 2.7 (a) if a < 1, then 

L(n) 

lim -— i-4 = 1. (2.10) 

n->oo log a - 



1 < nh(L(n)) and limsupn/i(L(?T,)) < — . (2-H) 

As in [13,14], we want to investigate a fluctuation process introducing a new "time" pa- 
rameter. In [13,14] this was done with a multiplicative parameter, but in our case an additive 
parameter is more adequate. This is caused by the hierarchical structure of the state space, 
which for a < 1 implies the logarithmic growth of the radius of the balls (Lemma 2.7(a)). Con- 
sider any non- decreasing function R : M + — > Z such that lim t _ ! . 0O R(t) = oo. For n = 1,2, .. . 
we define the fluctuation process as 

Y m N n ((L(n) + R(t))+) - EN n ((L(n) + R(t)) + ) f . 

X n\t) - 7== , l^-l^J 

(see q» . 

In what follows =^f denotes weak convergence of finite-dimensional distributions, and |_o;J 
is the integer part of x £ R. 

For a given by (12. 7p . we denote 

M-a 
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Theorem 2.8 Assume {2.7$ . 

(a) Suppose a < 1. Let (rij)j be any subsequence such that 

lim mh(L(ni)) = k (2.14) 

i— >oo 

for some n G [1,-]. Then X Ui =>f X^ K ' R \ where X^ K ' R > is a centered Gaussian process with 

EX^ R \s)X {K ' R) (t) = M R{sAt) g K (s V t), (2.15) 



covariance 



where 

_ [M-l) yj_ lye -e 

j=0 \k=0 / 

(b) If a = 1, iaen lim n _>. 0O nh(L(n)) = 1 and X n =>-y X^ R \ where X^ is a centered Gaussian 
process with covariance 

EXW(s)XW(t) = (1 - e- 2 )M R ( sAt l (2.17) 

Remark 2.9 (a) Existence of a subsequence (rij)j in part (a) follows immediately from Lemma 
2.7(b). For example, for c-rw we have L(n) = |l°g c /M and the condition (I2.14p is satisfied 
for any subsequence (n»)i such that log C//M ^- — |l°g c /M ^rj converges (to log c / M -), and any 
k 6 [1, — ) can be obtained in this way. 

(b) As it will be seen in the proof, part (b) of the theorem holds under the weaker (than 
rj+i/rj — > 1) assumption 

lim ^±1 = 0. (2.18) 

i^oo h[j) 

(c) It is well known that Gaussian processes with covariances of the form (12.151) have the Markov 
property. This is in sharp contrast with the corresponding result in the Euclidean case [13,14]. 

(d) The limit process obtained in part (b) has independent increments. 

(e) The function R(t) accounts for the time scaling of the limit process only. From the form of 
the covariances it is seen that the most natural forms of R(t) are 



R(t) 



ifa<l, (2.19) 



R(t) = \}og M t\ ifa = l. (2.20) 

We then have, for a < 1, 

t e < M R(t) < m e^ Q 21 j 

where 



9 = -^- T , (2.22) 



logM 
log ■ 
and 

t < M m <Mt if a = 1. (2.23) 
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It is also seen that the exponents in (I2.16P are "close" to ~, namely 

I < a m < 1. (2.24) 
t ~ at K J 

In (I2.2ip . (I2.23P and ( 12.24(1 . the left-hand inequalities are equalities for t = a~ n ,n e Z. 

Let us denote by and X the limit processes corresponding to R given by ( 12. 19(1 and 
(12.20(1 . respectively. Observe that in spite of the fact that such i?'s are not defined at t = 0, 
the processes X^ and X themselves are right continuous in L 2 at with X^(0) = and 
X(0) = 0. This follows from (12.211) . ( 12.23(1 . and since the function g K is bounded. 

Remark 2.10 The process X has the representation 



X t = VI - e- 2 W M n for M n <t< M n+1 , n G Z, 

where W is a standard Brownian motion. 

Properties of the process X^ are summarized in the next proposition. 

Proposition 2.11 (a) X^ is determined by a Gaussian sequence of random variables (Cn)nez, 

X { t K) = Cn fora- n <t<a-( n+1 \ 
where the ( n have the representation 



\n i / • i 



«=— oo 



with (i/i)igz iic?. standard normal. 

(b) X^ has a long range dependence property with exponent 1 m i/ie sense that 

< limsupr|E(X (K) (t) -X (K) (s))(X (K) (t + r) -X K (s + r))| < oo. (2.26) 

T — ^OO 

fcj For 6> > 1 (i.e., 1 > a > jj), if we define 

yW(t) = a m(e - 1)/2 xW(a- m t), t > 0, 

then 

a )lr a -n a -(n+i)) as m — ^ oo, 

where 

VeTT 



Y (t) = Z°±l f u e / 2 dW u , t > 0. 
t Jo 
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Remark 2.12 (a) Part (c) corresponds to Proposition 4.16 in [13]. Note that in that paper 
the limit process is a fractional Brownian motion, whereas in our case Y is a diffusion process 
satisfying the equation 

dY(t) = -~Y(t)dt + (V9~Ti)t e/2 - 1 dw t . 

Its covariance has the particularly simple form 

EY(s)Y(t) = (s A t) 9 (s V t)' 1 . (2.27) 

(b) The condition 9 > 1 in part (c) is required only to ensure that the process Y is well defined 
at 0. For 9 < 1 the assertion of part (c) remains true on [e, oo), e > 0. 

(c) We think that part (a) can be given a nicer form, analogous to that of the limit of Ym m 
part (c), i.e., that the process can be "interpolated" between the time points a n by a 
Gaussian diffusion process. To derive such a representation it would suffice to prove that the 
function g K defined by (I2.16p . with a R ^> replaced by ~, is decreasing in t. We have not been 
able to prove it. 

Remark 2.13 (a) Part (c) of Proposition 2.11 and Remark 2.12(b) show that the behavior 
of X^(t) as t — > oo changes drastically depending on whether 9 > 1 or 9 < 1 (for 9 > 
1, Hindoo VarX^^t) = oo, and for 9 < 1, lim^oo VaxX^it) = 0). In fact, the parameter 9 has 
a deeper probabilistic interpretation which is related to the recurrence/transience properties of 
the underlying random walk. These properties are in a sense characterized by the degree 7 of 
a random walk (see the Appendix). From (A. 2) it follows that 9 = 7 + 1. In particular, this 
implies that for 9 > 1 the random walk is transient, and for 9 < 1 it is recurrent. In the critical 
case 9 = 1 both situations can occur. Moreover, by Theorem 2.2 it is not hard to check that 
if 9 > 1, then lini^oo VarA^ n (L) = 00, and if 9 < 1, then lim^oo VarN n (L) = 0, while, as 
it was seen in Example 2.3, if 9 = 1 the number variance can be finite or infinite. We recall, 
however, that in general the condition lim sup^^ iV n (L) < 00 is not equivalent to recurrence 
of the random walk. 

(b) The equality 9 = 7 + 1 implies that for c-rw's the parameter 9 corresponds to - for the 
symmetric a-stable process in M (see (11. 3p . and (A. 2) in the Appendix). Hence for c-rw's the 
rate of the normine; ~ n 9 / 2 ) corresponds exactly to i 1 / 2 ", which is the norming in the 

fluctuation theorem in [13] (see (II. 5p ). 

Finally, we give the result for the Poisson system. 

Theorem 2.14 Under the assumptions of Proposition 2.5, for any rj-rw, any functions L : 
Z + — > Z + and R : K.+ — > Z, non- decreasing with lim n ^ 00 L{n) = 00 and lim^oo R(t) = 00, if 
X n is defined by K2~m) . then X n =^ CX^ R \ where X^ is given in Theorem 2.8(b) and C is 
a constant. 

Remark 2.15 From the proof of this theorem it will be seen that the same result is true if 
the initial number of particles at site u G Qm is T) u , where {jj u ) u ^ M are i.i.d. random variables 
with E(rj u ) 2+S < 00 for some 5 > 0, and Erf = Var r] u . 
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3 Proofs 

Let us denote 

B L (u) = {x G n M :\x-u\< L}, B L = B L (0), 

S L (u) = {x G ft M : \x - u\ = L}, S L = S L (0), 

p u n (L) = P( u + Z n eB L ). 

We write \A\ for the number of points in a bounded subset A of Qm (no confusion arises with 
the hierarchical distance). Note that 

\B L (u)\ = M L } \S L (u)\ = (M — l)M L -\ 

Observe that by ultrametricity of Qm, if u E B l , then u + £ n G B L if and only if £ n G B L . On 
the other hand, if \u\ = L + k, k — 1, 2, . . ., then B L C Si+^w), hence 

I -Br I 
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p*(L) = P(u + ( n eB L n S L+k (u)) = P(|e„| = L + k)- 

Thus we have 

| P(|£n| < L) if « g 5 L , 

Pl{L)= \ P ^\= L + k \ M -\)M^ if N = ^ + *,* = 1.2,.... (3 ' 1} 
Proof of Proposition 2.1 By fl2HD and M), 

EN n (L) = Pni L ) + E = 1^1" 

This proves (a) (this formula is surely known, but we have not found a reference ). 
In order to investigate VariV n (L) we study the tail of £„. We have 



P(\U>L) < Y / P(\£ 1 \<L,...,\Z k - 1 \<L,\£ k \>L) 

k=l 

n 

= Y, P t\p\<L) k - 1 P{\p\>L), (3.2) 

k=l 

by ultrametricity. Again by ultrametricity, the event that in the first n steps there is just one 
farthest jump whose length is j and it occurs at step k has the form 

A k,j = {|fi| < j, ■ ■ ■ , \€k-i\ < h = j, £fc+i e Bj^ k ), ...,£„ G Bj.x^fc)}, 

hence 



b(iu>^) > EE ^j) = E E < j) n ~ lp (\p\ = j) 

k=l j=L+l fc=l i=L+l 

> nP(|p| < L) n - x P(|p| > L). (3.3) 
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Formulas ( 13. 2 j) and (I3.3p imply that 

P(|£ n | > L) ~ > Ir) as L -> oo. (3.4) 

Using (12. ip and (13. ip we have 

VariV„(L) =/(L) + J(L), 

where 

J(L) = M L P(|£ n | < L)P(\U > L), 

m = m l p p(\u = L+k)(i- {M _ \ )Mk ^ P{\U = L+k)y 

Hence (b) follows immediately from ( 13.41) . 

The upper estimate in (c) is clear from the previous calculations. The lower bound is an 
easy consequence of 

VariV n (L) > K(L) > M L P(\£ n \ > L)(l - P(|&J > L)) 
and 1ET2]1 . ( 1331) . □ 

Proof of Theorem 2.2 All the statements follow from (I2.3P and the trivial formulas 

" r 1+L M^ L 



M L P(\p\ > L) = M L r L+1 + M L P(\p\ > L + 1), M L P(|p| > L) = ^ ^ fy 



□ 



Proof of Proposition 2.5 If the number of particles at each site is Poisson with parameter 
A, then 

VariV n (L) = \ PrX L ) = >^EN n (L) = XM L , 

by (12. 2p . hence the assertion follows. □ 

Proof of Lemma 2.7 (a) By ( 12. 7p and (12. 9 p it is easy to see that for any 0<e<aA(l — a) 
there exist positive constants Ci,C 2 such that for sufficiently large n, 

dia-ef <h(n) <C 2 {a + e) n . (3.5) 

Then, for large n, 



h 

which implies that 



1 , 1 

log a+e - + log a+£ — 

n c-2 



+ 1 ) < C 2 (a + e) log ^ 1/nC ' 2 = -, 

n 



L(n) < log a+e - + log a+e + 
n (_/2 
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hence 

L(n) log a 

lim sup y < 



Analogously, using the left-hand side of (I3.5P we obtain 

. log a 

hm mf r > 



n^oo l oga i log(a-e) 
Hence fl2TT0l follows 

(b) The lower bound in (12.111) is obvious by (12. 8ft . The assumption (12. 7p clearly implies 

lim l!i±* = a*, A; = 1,2,..., (3.6) 

n^oo r„ 

hence for a < 1 we have 

Y^(a + eV I < lim inf n+1 < lim sup n+ \ < ( i 

J n ^°° K n ) n^oo h{n) \j^ Q 

for any < e < 1 — a and k = 1, 2, . . .. The right-hand side inequality also holds for a = 1. 
Therefore, for all < a < 1, 

lim n+1 = 1 — a. (3.7) 

n->oo /i(n) 

By (J23]) and ([H]), 

nh{L(n)) < l + nr L(n)+1 , 
so f[2TTTj) follows easily from (JUD. □ 

Proof of Theorem 2.8 (a) We will prove 

lim Cov(X ni ( S ),X ni (t)) = M R ^g K {s V t) (3.8) 

and 

tLTM^wr- £ ^|^ )+ « (i))+ (- + o-^ + c<(^) + ^)) + ) 

(3.9) 

Then by (12.11) . (12 . 1 2[) and independence of random walks the result follows from the central 
limit theorem in the Lyapunov version. 
Denote 

L(i,t) = (L(ni) + R(t)) + . (3.10) 
Fix s < t. By (12. ip and independence of random walks, 

Cov(X ni (s),X ni (t)) = M~ L ^ ^p«(L(< >S ))(l-p«(L(i,t))) 

= M~ LM (I + E), (3.11) 



2+5 

= 0, <5>0. 
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where I = £ ueB{L(M)) • • • and J = E„ W( i, t )) • • •• We have 

/ = £ ...+ £ ... 

«eB(i(i,s)) u£B(L(i,t))\B(L(i,s)) 

= M L ^P(\U\ < L(t,s))P(\U > L(i,t)) 

+ E E p (i^=i)^f^(i-^i<L(M))) 



j=L(i,s)+l ueSj 

M L ^P(\^\ < L(i,t))P(\U > L(i,t)), (3.12) 



by (13. ip . Similarly 



' = E E 

j=i ue5 i(l>t)+j 



= M^) fjPd^l = +j) (l - F(|£J = L(i,t)+j) {M _ \ )W _^ j -(3.13) 

By (13.10n - fl3.13n . for sufficiently large i, we obtain 

Cov(X ni (s),X ni (t)) = M R ^\p(\U>L(htm + P(\U<L(i,t))) 



3=1 

It is known ([20], see also [7]) that 



1 f ^ 2 (IU = £(M)+j) l ( , u) 
M-l^ Mi- 1 J" K ' 



where 



Hence 



en 00 fn 

P{tn = u) = -jjfr + (M - 1) jfa> if |«| = *>0, 

j=k+l 



h = 1 - h(* - 1) - (3.15) 



\ 3=1 



M - 1 /j+fc+i /J+fc (3 16) 



M ^ Mi 

i=o 



This implies that 

M — 1 ,°° 1 

> L) = E ^(1 - /lVO- ( 3 - 17 ) 
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By ilTl . (TO , (BZZD and f!3~TUj) we have 

lim (h(L(i,t) +j - 1) + = Kte^" 1 , (3.18) 

4-S>00 y M — 1 / 

where 6 is defined by (I2.13p . Hence 
by <^m>- 

Combining (I3.16p - fl3.19p and f !3.14j) (it is clear that we can pass to the limits under the 
sums) we arrive at 

lim Cov(X ni (s),X ni (t)) 



M R{s) 



M — l °° 1 — e - KbaR<t)aJ \ ( M — 1 °° e - KbaR{t)aJ 



X} Mi 1 + M X} 



M z — ' M-? M M ^ Mi 

3=0 / V j=o 



(M — l) 3 °° 1 / °° e - Kbafl(t)ci:,+1+fc _ e -K.ba R ^ai X 2 



Mi 



M 4 ^ ikP' 1 ^ M k 

j=0 \k=0 

= M R( - s) g K (t) (3.20) 

(see (JUSD). This proves fETBl . 

It is easy to see that the expression under the limit on the left hand side of (13.91) is estimated 
from above by 

V V U ((L(n) +R(t)) + ) = - M {L{n)+R(t))+ 

M L(n)(l+5/2) Z-> PnKW') ^^)) ) M L(n)(l+5/2) 1V± 

by (Q, hence fl2~T01 follows. 

The proof of part (a) is complete, 
(b) The fact that 

lim n{h(L(n)) = 1 (3.21) 
follows immediately from (12. lip . This and (13. 7p imply that 

lim n (h((L(n) + R(t)) + + j - 1) + = 1 ^ 

n->oo y M — 1 / 

(cf. (BHD ), hence 

l l ^{L{n)+R{t)) + +j = e " 

for j = 0, 1, 2, . . .. Therefore, a counterpart of (13.201) is 

lim Cov(X n (s),X n (t)) = (1 - e^M^ 57 ^. (3.23) 

Now, the convergence (the Lyapunov condition (13.91) ) is obtained in the same way as before. 
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Observe that (I3.2ip . ( I3.22p . and hence (I3.23P as well, follow from (13. 7p in this case. (Recall 
that (13. 7p with a = 1 implies ( 12. lip , see the end of the proof of Lemma [2.7( b).) Therefore, as 
stated in Remark 2.9 (b), assumption (12. 7p can be replaced by ( I2.18p . □ 

Proof of Proposition 2.11 (a) From positive-definiteness of the covariance function ( I2.15P it 
is easy to see that M l / g K (a~ % ) is increasing in i, moreover, by (I2.16p . lim i ^_ 00 M l / g K (a~ l ) = 0. 
Hence the result is obtained by a direct computation. 

(b) By (I2.15p . fors<t<s + r<t + rwe have 

E(X {K \t) - X (K) (s))(X (K) (t + r) - X (K) (s + t)) 
= {M R ® - M R ^)(g K (t + r) - g K (s + r)), 

hence it suffices to investigate the second factor. Obviously, for large r, the time points s+r, t+r 
belong either to the same interval of the form [a~ k , a~( fc+1 )), or to two neighboring such intervals. 
Since, by (I2.19P and ( I2.16p . g K is constant on such intervals, it is enough to consider a sequence 
of the form 

T m = a~ m -d m , s < d m <t, m = 1, 2, 

Then, for large m we have 

T m (g K (T m + t) - g K {r m + s)) = T m a m a- m (g K (a m ) - ~g K (a m - 1 )), 

where g K (r) is obtained from ( I2.16P by putting r instead of a R ^\ The mean value theorem 
implies that \{g K {r) — g K ( L )) has a finite positive limit as r — )• 0. This proves (I2.26p . 

(c) Since the processes Ym^ are centered Gaussian, it suffices to prove convergence of covari- 
ances. By (I2TT9D and (T2~22D . 

R(cr m t) = m + R(t), m = l,2,..., (3.24) 

M m = a -m6 

This, (l2TT5j) and (T2~T9l imply that for s < t, 

Cov(yW( S ),yW(t)) = M R ^a- m g K (a- m t). 

By (I3.24p we also have 

_ e ~K,ba R ^' mt )ai 

lim = nba R ^a j . 

m->oo a m 

Hence, using the form of g K given by (I3.20p . and by ( 12 . 13j) . it is easy to see that 

lim Cov{Y m K \s), Y m K \t)) = 2KhI R{s) a R{t \ s<t. 

m— ¥oo 

Hence (c) follows (cf (|237jl ). □ 

Proof of Theorem 2.14 Let 7] u be the number of particles at site u at time 0. The r.v.'s 
(v u )uen M are i.i.d. Let (^' fc ) u gn A/ ,fceN be independent copies of r^-rw. Then 

N n( L ) = J2 J2 1 {u+tt' k eB L }' 

uGQm k=l 
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and for s < t, denoting L(n,t) = (L(n) + R(t)) + 

Cov(X n (s),X n (t)) = J2 [Ev u P(u + £ n eB L 



+ (Varrf - Er] u )P(u + ( n 6 B L{n , s) )P{u + £ n e B L(n , t) 
Er]° 



by Q. Hence 



lim Cov(X n (s),X n (t) = Er]°M R(s) . 

71— >00 

Next, using the inequality [a\ + . . . + a m ) 2+s < m 1+s (al +s + . . . + a 2 ^ 5 ), a, > 0, we easily obtain 

& IPM E E E 1b^, in + C k ) =0. 

«er2 fl/ \fe=i / 

Convergence of finite-dimensional distributions now follows from the central limit theorem 
(Lyapunov criterion). □ 

Appendix 



The development and applications of hierarchical random walks are outlined in [7] , and their 
recurrence/transience and some other properties are studied in [6J (and references therein). 

In order to deal with the standard a-stable process on W 1 and discrete time hierarchical 
random walks in a unified way, the continuous time version of a walk with unit rate holding 
time was taken in [5J. Its transition probability from to y in time t > is given by 

p t (0,y) = e-*y) -=?(»>((), y), 

n=0 

where p^(0,y) is the transition probability from to y in n steps. Generally, for a Levy 
process on a Polish space S (with additive group structure in the cases considered in [6j) with 
semigroup %, the degree of the process is defined by 

7 = sup{C > -1 : G c+ V < oo for all tp G B+(S)}, 

where B£(S) is the space of bounded non- negative measurable functions with bounded support 
on S, and 



i r°° 

G^ = Y(cjl t ^™ t ' c>0 ' 



is the fractional power of the potential (or Green) operator of the process. For £ = 7 it may 
happen that G 1+l ip < 00 or G 7+ V = oo,y? 7^ 0. Recurrence corresponds to 7 G (—1,0) and 
transience corresponds to 7 G (0, 00]. The value 7 = is special since both cases G<p < 00 and 
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G(p = oo, ip 0, can happen. 7 = 00 occurs, for example, for a simple asymmetric random 
walk on Z. For 7 > 0, 7 is also given by 



7 = sup{C > : EL Br < 00 for all R > 0}, 

where Lb r is the last exit time of the process (starting at 0) from an open ball Br of radius R 
centered at 0. For a discrete time random walk on a discrete space, instead of G^ip it is natural 
to consider 

_£y-y»>(o,o), (a.i) 

^(o,o)^=f; r(n+ f + V >(o,o), 



but since 



n! 

n=0 



and 

r(n + c + i) , 

— >■ 1 as n — )■ 00, 

nivS 

the degree of the random walk is the same in discrete and continuous time. 

For the standard a-stable process on M. d , 7 is given by (11. 3p . and the range of possible 
values of 7 for those processes is restricted to the interval [— |, 00). For the c-rw defined 
by (II. II) . (11. 2p . 7 is given by (11.41) . and the range of possible values of 7 is (— l,oo). Hence 
this class of hierarchical random walks is richer than the class of standard a-stable processes 
and corresponding symmetric random walks on Z d . The c-rw's mimic the behavior of a-stable 
processes in the sense that they have the same recurrence/transience properties for equal values 
of their degrees; in particular, for 7 = 0, dp = 00 holds for both of them. The correspondence 
of degrees allows to choose c in order to study "caricatures" of a-stable processes by means of 
c-rw's, including non-integer values of the dimension d (which is one of the reasons for using 
hierarchical random walks in statistical physics, see [8] and references therein). But there 
are also differences, for example, the distance from of the c-rw has a different behavior from 
the a-stable Bessel process (see [6], Remark 3.5.7). The correspondence between the degrees of 
the a-stable process and the c-rw also serves to compare our results with those of the a-stable 
process and exhibit analogies and differences. 

The assumption rj +1 /rj — >aasj— >oo,0<a< 1, implies that the degree of the r 3 -rw is 
given by 

logM 

7 = 5 1 for a < 1, 7 = 00 for a = 1. (A. 2) 

log - 

This is obtained from [6] (Proposition 3.2.7). Hence 

< -jjjr, recurrent 
j_ 

M ' 

> jj, transient. 

The parameter 9 = 7 + I (see (I2.22p ) therefore has a meaning in terms of recurrence/transience. 
Examples: 
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1. c-rw: < c < M, 

•H 1 -£)(£)'"• ' = 1 - 2 



c \ / c y"-i c logc 
— II — I 1 = 1.2,..., a = — , 7 = 



M' ' bg(M/c)' 
recurrent for c < 1, transient for c > 1. 
2. /-nu:/3 > 0, 



where D is a normalizing constant, 



a = W 7 = ' 



recurrent for /3 < 1, transient for /3 > 1 (follows from the recurrence criterion in [TT], see 
also Example 3.2.6 in [6], which is slightly different) . 

3. 

r j = Dj-^ j = l,2,..., 
where > 1 and Z) is a normalizing constant, 

a = 1, 7 = oo, transient. 



Remark. In |6j rj is written in the form 

rj = jD ivfe' ^' = 1 ' 2 >---> 

where is a positive constant, (cj) is a sequence of positive numbers, and D is a normalizing 
constant. The parameter fi is useful for some applications (see also [3 IB])- 
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